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In a previous paper, we introduced a new interpretation of matrix models, in which any 
d-dimensional curved space can be realized in terms of d matrices, and the diffeomorphism 
and the local Lorentz symmetries are included in the ordinary unitary symmetry of the 
matrix model. Furthermore, we showed that the Einstein equation is naturally obtained, if 
we employ the standard form of the action, S = —tr ([A a: At][A a , A b fj + ■ ■ ■ . In this paper, we 
extend this formalism to include supergravity. We show that the supercovariant derivatives 
on any d-dimensional curved space can be expressed in terms of d supermatrices, and the 
local supersymmetry can be regarded as a part of the superunitary symmetry. We further 
show that the Einstein and Rarita-Schwinger equations are compatible with the supermatrix 
generalization of the standard action. 



■ 81. Introduction 

o 



Although it is believed that string theory can provide a formalism describing 
the unification of the fundamental interactions, its present formulation based on 
O , perturbation theory is not satisfactory. In order to examine whether it actually 

describes our four-dimensional world, a non-perturbative, background independent 
formulation is needed. Matrix models represent a promising approach to the study 
of the nonperturbative dynamics of string theory. For critical strings, they are ba- 
sically obtained through dimensional reduction of the ten-dimensional U(N) N = 1 
supersymmetric Yang-Mills theory. 2 )' 3 ) 

IIB matrix model 3 ) is obtained through dimensional reduction to a point, and 
the action is given by 

S = ~tr (j[A a ,A b ][A a ,A b ) + ^ 7 a [A a ,^]\ , (1-1) 

where ip is a ten-dimensional Majorana-Weyl spinor, and A a and ip are N x N 
Hermitian matrices. The indices o and b are contracted by the flat metric. This 
action has an 5*0(10) global Lorentz symmetry and U(N) symmetry. A problem 
with this model is that it is unclear how curved spaces are described and how the 
fundamental principle of general relativity is realized in it. 

It may seem that the space of dynamical variables would become very small 
through the dimensional reduction, but this is not the case if N is infinitely large. 4 ) 
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Indeed, IIB matrix model contains super Yang-Mills theory: 

dim red N — >oo 

10D super Yang-Mills — -> ' IIB matrix model D 10D super Yang- Mills. 

The matrix-valued variables A a and V'a act on the Hilbert space V = C N as endo- 
morphisms, i.e. linear maps from V to itself. Because V is infinite dimensional in the 
large- iV limit, we can interpret it in various ways. If we assume that V is the space 
of an n-component complex scalar field, V = {<£ l : M 10 — > C n }, instead of V = C N , 
then an endomorphism T is a bilocal field K lJ (x, y), which can be formally regarded 
as a set of differential operators of arbitrary rank with n x n matrix coefficients: 5 ^ 

(Ti P ) i (x)=Y i / d 10 yK^(x,y)^(y) 

= cg } (x)<f? (x) + cff (x)d^ (x) + c^j ij (x)d^ (x) + • • • . (1-2) 

In particular, we can consider the covariant derivative as a special value of A a : 

A a = i(d a - ia a (x)) G End{V). (1-3) 

In this sense, the ten-dimensional U(n) M = 1 super Yang-Mills theory can be 
embedded in the matrix model, and local gauge symmetry is realized as a part of 
the U(N) symmetry of the matrix model; that is, we have 

5A a = i[X,A a ], (1-4) 

where A is a matrix- valued function on M 10 , which is a 0-th order differential operator 
in End{V). This interpretation of V, however, does not manifest the existence of 
gravity. Therefore it is desirable to formulate another interpretation in which the 
diffeomorphism and the local Lorentz symmetries are embedded in U(N). With such 
an interpretation, any curved space corresponds to a certain matrix configuration, 
and the path integral includes the summation of all the curved spaces. Our idea is 
that it is not a curved space itself but, rather, the covariant derivatives on it that 
correspond to matrices. Note that all the information needed to describe the physics 
on a manifold are contained in the covariant derivatives. 

Suppose we have a curved space M with a fixed spin structure and a covariant 
derivative on it*' . Our task is to find 

(a) a good space V and 

(b) a good object V( ) that is equivalent to a covariant derivative V a 

such that each component of V( a ) (a = 1, 2, • • • , 10) is expressed as an endomorphism 
on V. In a previous paper, 1 ) we showed that V is given by the space of functions on 
the principal Spin(10) bundle over M and that V( tt ) is given by 

V (o ) = R(a)\g~ 1 )^b, (1-5) 



*' Because a covariant derivative introduces a new vector index, it is not an endomorphism. 
Therefore, it cannot be represented by a set of matrices. 
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where V a is the covariant derivative 

V a = e a m (x) (d m - u m bc (x)O bc ) , (1-6) 

and R( a } b (g) is the vector representation of Spin(10). Here we assume that all indices 
are Lorentz. We can show that each component of V( a ) is indeed an endomorphism 
on V, while V a is not. In we review this formulation in detail. By introducing 
the space V given above, we can successfully embed gravity in the matrix model in 
such a way that the diffeomorphism and the local Lorentz symmetry are realized as 
parts of the U(N) symmetry of the matrix model. 

In contrast to the situation with gravity, however, supergravity cannot be easily 
embedded in the matrix model in this interpretation, because the U(N) symmetry 
does not contain a fermionic symmetry. In order to implement a manifest local 
supersymmetry in the matrix model, we need to extend the space V and the U(N) 
symmetry to include a fermionic parameter. This is done by extending M to a curved 
superspace and taking V to be the space of functions on the principal Spin(W) 
bundle over the superspace. This extension forces us to consider supermatrices 
instead of ordinary matrices. 

The organization of this paper is as follows. In the next section, we review our 
previous paper. 1 ) First we explain how a covariant derivative on a d-dimensional 
curved space can be expressed as a set of d endomorphisms. Then we introduce a 
new interpretation in which matrices represent such differential operators. Based on 
this interpretation, we can show that the Einstein equation follows from the equation 
of motion of the matrix model. In we express a supercovariant derivative on a 
curved superspace in terms of endomorphisms. Although we consider only TV = 1 
supersymmetry, the generalization to Af > 1 is straightforward. In we introduce 
supermatrix models. Dynamical variables are regarded as supercovariant derivatives 
on a curved superspace, and the symmetries of supergravity are embedded into the 
superunitary symmetry. We show that the equations of motion of the supermatrix 
model are satisfied if we assume the standard torsion constraints and the supergravity 
equations of motion. JDis devoted to conclusions and discussion. In Appendix IA1 we 
summarize the properties of supermatrices. In Appendices [B] and [O we present the 
detailed calculations in deriving a classical solution of the supermatrix model. 

§2. Describing curved spaces by matrices 

In this section, we summarize the results given in Ref. 1). We first explain how 
a covariant derivative on a cf-dimensional manifold M can be expressed by a set of d 
endomorphisms acting on the space of functions on the principal Spin(d)-bimd\e over 
M. We then apply this idea to the matrix model and introduce a new interpretation 
in which dynamical variables represent differential operators on curved spaces. We 
show that the matrix model reproduces the Einstein equation correctly, and that the 
symmetries of general relativity are realized as parts of the U(N) symmetry. 
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2.1. Covariant derivative as a set of endomorphisms 

Let M be a Riemannian manifold with a fixed spin structure and M = UiUi be 
its open covering. On each patch U, the covariant derivative is expressed as 

VH = e a m ^(x) (d m - u m bc ^(x)O bc ) , (2-1) 

where e a m (x) and u> m bc (x) are the vielbein and the spin connection, respectively. O bc 
is the Lorentz generator that acts on Lorentz indices. The index [i] is the label of 
the patch. In the overlapping region U n Uj, the operators and are related 
by 

V^=R a b ( tlJ (x))V b b \ (2-2) 

where tij : U{ D Uj — s- G = Spin(d) is the transition function, and R a b (tij(x)) is the 
vector representation of tij(x). 

Let us consider the principal G bundle on M associated with the spin structure 
and denote it by E pr i n . It is constructed from the set Ui x G by identifying (xm, gui) 
with (x m ,g y] ): 

x [i\= x \j], 9[i\=Uj{x)g\j]- (2-3) 

We take V = C 00 (E pr i n ), which is the space of smooth functions on E pr i n . We 
assume that covariant derivatives act on the space V, that is, O ab generates an 
infinitesimal left action, 

ie ab (o o6 / [i] ) (x, g) = U(l + ie ab M ab ) ~' g\ - /M (x, g) , (24) 

where M ab is the matrix of the fundamental representation. Then, we can construct 
endomorphisms from a covariant derivative in the following way: 

vfa) = *k*)%f)vf- ( 2 - 5 ) 

Here, R^ b (g) is the vector representation of G.*^ In the overlapping region of two 
patches, the actions of and V/jl on / G V = C 00 (E pr i n ) are related as 

vjy^M,^) = R {a) %i)V b i] f®(x[H,9\i]) 

= R [a) b {g^)R h c {t l] {x))vffi\x m ,g [3] ) 
= R [a) c ((^(x)- 1 ^)- 1 ) VffW{x m ,g m ) 
= R {a) c (g b l)v^f l Hx m ,g {j] ) 
=vg ) / w (x w> p M ) ) (2-6) 

*' R a b and -fya) 6 represent the same quantity. However, we distinguish them, because a and (a) 
obey different transformation laws. Specifically, a is transformed by the action of G, while (a) is 
not. 
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where f^ % \xuu gu\) is the expression of / on Ui x G. Note that with the identification 
given in ()2-3|) . we have /^(sm, flta) = f^( x \j]i 9]j])i by definition. This confirms that 
each component of V( a ) is globally defined on E pr i n and is indeed an endomorphism 
on V. The index (a) merely labels d endomorphisms. Similarly, operators with 
vector or spinor indices can be mapped to a set of endomorphisms using the vector 
and spinor representation matrices of G, for example, 

T ia) = R (a) p (g- l )T^ 
T (a){b) = R(a) c {9~ 1 )R(b) d (9~ 1 )T c d- (2-7) 

Here i?( a ) b (g _1 ) and R^ /3 (g~ l ) are the vector and spinor representations, respec- 
tively. In addition, we have the relation 

V ia) f b (x,g) = R b ic) (g)V {a) f {c) (x,g), (2-8) 

where the covariant derivative on the left-hand side also acts on the vector index b. 
This point is explained in detail in Ref. 1). 

The method presented in this subsection is valid in any number of dimensions, 
and we can express the covariant derivative on any ci-dimensional Riemannian man- 
ifold in terms of d matrices. 

2.2. A new interpretation of IIB matrix model 

In this subsection we present a new interpretation of IIB matrix model. ^ As 
we showed in the previous subsection, for any <i-dimensional Riemannian manifold 
with a fixed spin structure, its covariant derivative can be described by a set of d 
matrices acting on C 00 {E pr i rt ). 

Let us consider the space of large- N matrices. For any manifold M, this space 
contains a set of matrices which are unitary equivalent to the covariant derivatives 
iV( a ) on this manifold. If the matrices A a are sufficiently close to the derivatives 
iV( a ) on one of the manifolds M, it is natural to regard these A a as acting on 
C°° (Ep r i n ), and to expand A a about iV( ):*) 

A a =*V( a ) + a {a) (x,g) +ia (a) (6) (x,2)V (6) + ia {a) bc (x, g)O hc 

+i 2 a (a) ^ c Hx,g)V {b) V {c) +i 2 a {a) ^ cd (x,g)V {b) O cd + -- - . (2-10) 

In this expansion, local fields appear as coefficients. For example, a^ b \x, g) con- 
tains a fluctuation of the vielbein. Coefficients of higher-order derivative terms cor- 
respond to higher-spin fields. In this sense, this part of the space of large- N matrices 
describes the dynamics around the background spacetime M. Some matrix configu- 
rations correspond to S* 10 , others to T 10 , and so on (see Fig^)**). 

*' Strictly speaking, because A a is Hermitian, we should introduce the anticommutator { , } in 
Eq. (HUSt : 

A a = iV( o) +a (a) (x,g) + l -{a (a) w (x,g), V (6) } + ^{a (a) bc (x, g), O bc } H . (2-9) 

In the following, we omit the anticommutator for simplicity. 

**' A matrix configuration in which some of matrices A a fluctuate around corresponds to a 
lower-dimensional manifold. 
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Note that we do not fix M and that all Riemannian manifolds with all possible 
spin structures are included in the path integral. Of course, it also contains matrix 
configurations that are not related to any manifold. For example, the matrices A a 
satisfying [A a , Af,] = iB a i,-l with constant B a f, describe a flat noncommutative space. 
Another trivial example is A a = (o = 1, 2, • • • , 10). 



Classical solutions 

Now let us consider the equation of motion. The action of IIB matrix model is 
given by Eq. (|1T(I . Variation with respect to A a gives the equation of motion 



Not manifolds 



Space of large N matrices 




Fig. 1. Path integral include all the curved spaces. 



[A a ,[A a ,A b }] = 0. 



(2-11) 



Here, we simply set ip = 0. If we impose the ansatz 



K = iV( a ) 



(2-12) 



then Eq. (|2-ll|l becomes 



v (a) ,[v (a) ,v (6) ] =0. 



(2-13) 



Equation (|2-13l) is equivalent to*) 



= [V 
= [V 



a 



a 



[V Q ,V 6 ]] 

Rab Cd O c d\ 




(2-14) 



We use the formula 12-81 and T R = R 



-i 
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where R a b cd is the Riemann tensor and R b c = R a b ac is the Ricci tensor. Here we 
have assumed that V a is torsionless. Note that Eq. (|2- 14f) holds if and only if we 
have 

V a R ab cd = 0, R ab = 0. (2-15) 

The first equation here follows from the second by the Bianchi identity V^ a R bc de ^ = 0. 
Therefore, the covariant derivative on a Ricci-flat spacetime is a classical solution. 
We next consider a matrix model with a mass term: 

5 = ~]f tr (j[Aa,A b ][A a ,A b ] + i^y°[4»,^A + ^tr(A a A a ) . (2-16) 

The equation of motion is now given by 

[A b ,[A a ,A b ]]-2m 2 A a = 0, (2-17) 

which leads to 

R a b = -2m 2 5 a b . (2-18) 

This is the Einstein equation with a cosmological constant 

A = -8m 2 . (2-19) 

However, if m 2 > 0, the configuration A a = (a = 1, ■ • • , 10) is also a classical 
solution and has lower free energy. Therefore, such a configuration may dominate 
the path integral. If the positive mass term exists only for d directions, 5 mass = 
^p- Yla=i ^ r {A a A a ), then we expect that these d directions are collapsed to a point 
and the rest 10 — d directions remain flat. 

The analysis of the eigenvalue distribution 6 ) is consistent with this expectation. 
Suppose an effective mass term is generated in six directions and not in the rest 
four directions. Then that analysis indicates that the four dimensional spacetime is 
generated. 

DifFeomorphism and local Lorentz invariance 

We now show that the symmetries of general relativity are realized as parts of 
the U (N) symmetry of the matrix model. The U(N) symmetry of the matrix model 
is expressed as 

8A a = i[A,A a ], (2-20) 
where A is a Hermitian matrix. If we expand A a as Eq. (j2- lUf) and take 

A= l -{X^\V {a) }= l -{X a (x),V a }, (2-21) 

Eq. (|2-2U|) becomes a diffeomorphism on the fields that appear in Eq. (|2-10|) . and 

e a m and uj m ab that appear in V( a ). For example, we have 

5e a m (x) = e a n (x)V n X m (x), (2-22) 
5oj m ab {x) = X n (x)R mn ab (x), (2-23) 
5a a (x) = -X n V n a a {x). (2-24) 
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This is the transformation law of the diffeomorphism. Similarly, if we take 

A = i\ ab (x)O ab , (2-25) 

we obtain 

5e a m {x) = -X a \x)e b m {x), (2-26) 
Su m ab (x) = -d m X ab (x) + 2\t a c (x)u Jm b ^(x), (2-27) 

which is the local Lorentz transformation. 

§3. Supercovariant derivative as a set of endomorphisms 

In the last section we showed that the symmetries of general relativity can be 
embedded in the U(N) symmetry of the matrix model. In order to implement local 
supersymmetry, we employ a superspace M as the base manifold and consider the 
supercovariant derivative. 

3.1. General prescription 

In this subsection, we introduce the supercovariant derivative on a d-dimensional 
M = 1 superspace and explain how to relate it to a set of endomorphisms. 

The coordinates of the superspace are z M = (x m , 9^), where x m (m = 1, 2, • • • , d) 
are the bosonic components and 9^ {/jl = 1, 2, • ■ ■ , d s ) are the fermionic components. 
The quantity d s is the dimension of the spinor representation. For example, if d = 10, 
9 corresponds to the Majorana-Weyl spinor and d s = 16. If d = 4, 9 corresponds to 
the Majorana spinor and d s = 4. Also, M = (m, /i) is an Einstein index, which can 
be converted into a Lorentz index, A = (a, a), using the supervielbein E^ M (x, 9). In 
the following, we denote the Einstein indices by M, N,L, - ■ ■ and the Lorentz indices 
by A,B,C,---. 

The supercovariant derivative V which acts on Lorentz indices is defined as 

V M = 9m - ^M^Oab, (3-1) 

where dyi = (g^m, gijr) is the left-derivative and O ab is the Lorentz generator. Note 
that O a b has only Lorentz-vector indices. The action of V on an Einstein index is 
defined as follows: 

V M V N = (-) m ^E N B V M VB 

= d M v N + r MN L v L . (3-2) 

Here, m, n and b are if M, N and B are bosonic indices and 1 if they are fermionic 
indices. The quantity Fmn L is the Christoffel symbol, which is chosen to satisfy 

V M E N A = 0. (3-3) 

The commutation relation of = Ea Vm is given by 



[V A , V b } = R AB cd O cd + T AB C V C , 



(3-4) 
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where Rab is the super Riemann tensor and Tab is the supertorsion. 

In H2.ll it is shown that the covariant derivative on a curved space can be 
expressed by a set of matrices. In the same way, the supercovariant derivative on a 
superspace can be expressed by a set of supermatrices. 

Consider the principal G = Spin(d) bundle £ pr i n over the superspace M.. If M. = 
UiUi is an open covering and Z[i] = (x^,0^) and zy] = (xy],6\j]) are the coordinates 
of the same point in UiDUj, then £ pr in is constructed from Hi x Spin(d) by identifying 
( z [i\,9[i\) an d (%])%]) i n the case ^ = Uj(z)gy], where %(z) is the transition 
function. Taking V = C°°(£p r j n ) = {/ : £ p ri n — > C}, we can convert the Lorentz 
indices a and a of the supercovariant derivatives into those with parentheses :*) 

V( o) = R {a) b (g^)V b , V (a) = R^fig-^Vp. (3-5) 

As in the derivatives V( a ) and V( a ) are globally defined on £ pr i n and, indeed, 
are endomorphisms on C°°(£ pr i n ). The only difference between the two cases is that 
in the present case, transition functions, parameters of diffeomorphisms and local 
Lorentz transformations depend not only on the bosonic coordinates x m but also on 
the fermionic coordinates 9^. More specifically, V( a ) acts on / £ V = C°°{£p r i n ) as 
follows: 

v[VM(^ ] , 5[i] ) = i? (a) 6 ( 5[ T ] 1 )v[ i] / [il (^],5 M ) 

= i? (o) 6 ( 5 - 1 )i? 6 c (^(z))vW/ b1 (%],^l) 
= R {a) c m j ( z )- 1 g [{l y 1 )V^f^(zy ] ,gy ] ) 
= R {a f(g^)vWfW( Zlj] ,g m ) 

= ^} ) f lj] (^m)- (3-6) 

Similarly, we also have 

v[V H (^]^ H ) = vf2 ) /W( %] , 5b1 ). (3-7) 

Therefore, the V( a ) and V( Q ) are endomorphisms on a supervector space C°°(£p r i n ), 
which are described by supermatrices, as we explain in the next section. 

Before concluding this subsection, we give a proof that V( ) and V( Q ) are anti- 
Hermitian and Hermitian, respectively. The standard inner product is defined by 

(u{z,g),v(z,g)) = J dg J Ed d xd d °6u*(z,g)v(z,g), (3-8) 

where dg is the Haar measure of Spin(d) and E = sdetEu A - Then, we can show 
the (anti-) Hermiticity of V^i in the following way: 

(u [a) {z,g)M a) v{z,g)) = j Ed d xd d °6ul a) (z,g)V^v(z,g) 

*) jj(spinor) /3 - g ^ ne a pp r0 p r i a t e spinor representation. It is Majorana-Weyl for d = 10 and 
Majorana for d = 4. 
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Ed d xd ds 9u* a (z,g)V a v{z,g) 

Ed d xd d °9(V a u* a (z,g))v(z,g) 

= - f Ed d xd d °9 (vW^ a) (*,<?)) v(z,g) 

= -(v^u {a) (z,g),v(z,g)y (3-9) 

Similarly, we can show that 

(u (a) (z,g)M a) v(z,g)) = (v^u {a) (z, g),v(z, <?)) . (3-10) 

The extra minus sign appears because u a is fermionic. 

The method presented in this subsection can be applied to a supersymmetry 
with any value of M. 

3.2. Example:two- dimensional Af = (1, 1) superspace 

As a concrete example, let us construct a two-dimensional M = (1,1) superspace 
M with Lorentzian signature and a principal Spin(l, 1) bundle £ pr i n over it. 

In order to eliminate the fields that are unnecessary in supergravity, we impose 
the following torsion constraints: 7 ) 

T ab c = T af ? = 0, T a(3 c = 2i(C-f c ) af3 . (3-11) 

Here C a g is the charge-conjugation matrix. Let M be the ordinary manifold which 
is obtained by setting all the fermionic coordinates 9 of M to zero. We denote by 
e a m (x) and u m (x) = ui m m (x) the zweibein on M and the torsionless spin connection 
associated with e a m (x), respectively. Then, using the higher-^ components of the 
superdiffeomorphism symmetry, we can write the superzweibein E(x, 9) and the spin 
connection f2(x,9) in the form of the Wess-Zumino gauge: 7 ) 

pi A i m _ ( e m a (x) -\u m {x){ioi0) a 

Em {x ' 9) -{ i (c 1 «e), 
F M (rfj)= ( e a m {x) \u a { x ){ lm ey \ , , 

n m 01 (x,e) = n m (x,d)=u Jm (x), nf{x,o) = n lx (x,e) = n. (3-13) 

Here we have set the gravitino and auxiliary fields to zero for simplicity. The super- 
covariant derivative Va = Ea m Vm is given by 

V a = E a M X7 M = e a m d m - 2u a \o 01 - 7(70^)^) , (3-14) 



4 



V a = E a M V M = Scfdp - i(C 7 b e) a V b . (3-15) 

Although Oqi does not act on the index /1, if we regard C01 — 3(701^)^^ as a "Lorentz 
generator" , it does act on fj,, and therefore 9^ can be regarded as a "Lorentz spinor" . 
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Then V a can be regarded as the ordinary covariant derivative on M which also acts 
on the index fi. In this sense, A4 is a spinor bundle over M with Grassmann-odd 
fiber coordinates 6^, which is associated with the spin structure. 

Now we introduce a principal Spin(l, 1) bundle £ pr i n over A4 associated with 
the spin structure. Then, we regard the supercovariant derivative as acting on 
£prin- Parameterizing the spin s representation of Spin(l, 1) as 

R( s )(<p) = e 2sv , (3-16) 

we can express the Lorentz generator as 

Ooi = \d v . (3-17) 

Then, if we introduce 0^' according to the relation 

qO*) = ( e ^7oi)(M) i/ ^ ; ( 3 .18) 

qM i s globally defined on £ pr i n . Next, using x m , eM and tp as independent variables, 
we can rewrite Eqs. (|3-14(l and 1)3-15(1 as 

V a = e a m (d m - ^Wmdtp) , 

V a = (e^)^ a d M - i(C 7 fe e^) Q V . (3-19) 

Finally, as explained in the previous subsection, we can convert indices without 
parentheses into those with parentheses, and hence we have 

V (o) = e {a) m (d m - ^uJ m dA , (3-20) 



V (a) = *(a) (M) d M - i ( C 7 (6) <?) V (6) . (3-21) 



(a) 

Note that £ pr %n can be constructed direct product of M°l 2 = {8^,9^} and E pHn . 
Supersphere 

We now construct the homogeneous and isotropic Euclidean supersphere S 2 ' 2 
and the covariant derivative on it. Let us start with the case of the ordinary S* 2 . 1 ) 
Then, by parameterizing Spin(2) as 

R( s \ l p)=e 2isv , [0,2vr)) (3-22) 

we can express the Lorentz generator in terms of the derivative with respect to ip as 

0+- = (3-23) 

where + and — indicate the linear combinations of the Lorentz indices 1 + i2 and 
1 — i2, respectively. 
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In the stereographic coordinates (w, w) projected from the north pole, the metric 

is 

9ww — 9ww — ~jz ', -\o i 9ww — 9ww — 0, (3 - 24) 
(1 + ww) z 

and we can take the zweibein and spin connection as 

1 

e +w = e-w = 0, (3-25) 



1 + WW 



u w = z ' -, = -z—. -■ (3-26) 

1 + WW 1 + WW 

Therefore, the covariant derivatives with parentheses are given by 

= e" 2 ^ ((1 + ww)d w + l -wd^j , 

= e 2i ^ ((1 + ww)d iij - l -wd^j . (3-27) 

In the same way, in the stereographic coordinates (w',w') projected from the south 
pole, which are related to (w,w) as w' = 1/w, we have 

V$ = e" 2 ^' ((1 + w'w>)d w , + i*fdj\ , 

VK ] = e 2 ^' ( (1 + w'w')d^ - l -w'd v ) . (3-28) 



(_) — c ~ w w )"v)' 2 

Here, the coordinate tp' of the fiber on this patch is related to <p according to the 
transition function 

^ = V9 + argH + |. (3-29) 

We can explicitly check that Eqs. ()3-27[) and H3-28[) are identical. Note that E pr i n 
on S" 2 is Spin(2)(= S 1 ) bundle over S 2 and is topologically S 3 . 

Next we consider the supersphere S* 2 ' 2 and the principal Spin(2) bundle £ pr in 
over it. In this case Eqs. (|3-2Uj) and ()3-21|) become 

V M = | (1 + w - )Qw + * 

V M = e 2i V | (1 + W - )Q _ _ i_ wd \ f (3 . 30) 
v (5=*(«) M »M- i ( G T' (6)fl ) (tt) v W- ( 3 - 31 ) 



Using the relation (w,9ri) = (l/w',9^) and Eq. (l3~m we can explicitly confirm 
the relation 



V^ = Vg. (3-32) 
In this way, we can describe the supersphere 5 2 ' 2 in terms of four supermatrices. 
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§4. Supermatrix model and local supersymmetry 

In the previous section, we showed that the supercovariant derivative on a su- 
perspace can be described by a set of supermatrices. Therefore, as in ^2.21 we can 
embed local supersymmetry manifestly in supermatrix models. In this section, we 
consider a straightforward generalization of IIB matrix model. We show that the 
equation of motion of the supermatrix model is compatible with that of supergravity. 

4.1. Supermatrix model 

An element / : £ pr i n — > C of V = C°°{£pri n ) can be expressed as a power series in 
0L l . The coefficients of even powers of 6^ in this series are Grassmann-even functions 
and form the Grassmann-even subspace V e of V. Similarly, the coefficients of odd 
powers of 6^ form the Grassmann-odd subspace V Q . Therefore, the total space V is 
the direct sum of V e and V a : 

V = V e ® V . (4-1) 

If we introduce a regularization of V in which the dimensions of V e and V are 
N e and N a , respectively, then a linear transformation on V can be expressed as a 
(N e , A r Q )-supermatrix. We summarize the definition and properties of supermatrix 
in Appendix !a1 

The derivatives V7 a ) and O a b map an even (resp., odd) element to an even (resp., 
odd) element. Hence they can be expressed in terms of even supermatrices. Con- 
trastingly, V( a ) maps an even element to an odd element and vice versa. Therefore 
they are expressed in terms of odd supermatrices. 

Now we propose the supermatrix generalization of IIB matrix model. The dy- 
namical variables consist of Hermitian even supermatrices A a with vector indices 
and the Hermitian odd supermatrices W a with spinor indices. The action is given by 

S = ~j 2 Str (^[A a ,A b )[A\A h ) + \* a {^JiAa^ . (4-2) 

This model possesses superunitary symmetry U(N e \N ). If the matrices A a are 
sufficiently close to the covariant derivatives «V( a ) on some £prim we regard the A a 
and ty a as endomorphisms of V = C°°(£p r i n ). As we discuss in the next section, 
there is a possibility that Eq. (|4-2j) is equivalent to the original IIB matrix model. 

A remark on the role of the odd supermatrix & a is in order here. The action 
(|4-2|) has the global symmetry 

6^Aa = -i^lae, 5^ a = l -[A a ,A b ][r b e) a (4-3) 

and 

5^A a = 0, 6®V a = e«, (44) 

which is the analogue of the J\f = 2 global supersymmetry of the original IIB matrix 
model. Here, the quantities 
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are odd scalar supermatrices, which commute with A a and anticommute with & a . 
In the original interpretation of IIB matrix model, this symmetry is simply the ten- 
dimensional N = 2 supersymmetry. In the case of the supermatrix model, however, 
the local supersymmetry is manifestly embedded as a part of the superunitary sym- 
metry, as we see in N4.3I Indeed, a supermatrix model without such as 

S = ~^Str ([A a ,A b ][A a ,A b ]) , (4-6) 

also possesses this local supersymmetry. Furthermore, as we see in the next sub- 
section, both Eqs. ()4-2[) and (|4-6p describe supergravity at the level of classical 
solutions. On the other hand, in the original interpretation of IIB matrix model, 
in which matrices are regarded as coordinates, the TV = 2 global supersymmetry 
prevents the eigenvalues from collapsing to one point, and this allows a reasonable 
interpretation of spacetime. 3 ) We conjecture that the global symmetry represented 
by Eqs. (|4-3j) and ()4-4|) plays a similar role here and allows A a to fluctuate around 
*V( ) on some manifold. 

4.2. Classical solutions 

In this subsection we investigate classical solutions of Eq. (|4-2|) . For simplicity, 
let us consider the four-dimensional model here. We use the notation of Ref. 8). In 
this case, A a and \P a are four-dimensional vector and Majorana spinor, respectively. 
The equations of motion are given by 

[A a , [A a , A b \] + {C lb ) a P{* a , 9 P } = 0, (4-7) 
(7 a )/[A,^] = 0, (4-8) 

where C a @ is the charge-conjugation matrix. 
We impose the following ansatz:*) 

A a = iV {a) , <F Q = 0. (4-9) 

In the bosonic case, once we impose the ansatz (|2-12|) . we can obtain the general 
solution of the matrix equation of motion ()2-llj) . In the present case, however, it 
is not easy to find the general solution of Eqs. ()4-7|) and (|4-8(l . Instead, we show 
that the ansatz (|4-9|) satisfies Eqs. (|4-7|) and (|4-8|) if we impose the standard off-shell 
torsion constraints 9 ^' 10 ^ 

V = 2i(7 c C~ 1 ) Q/3 , 

If we start with A a = iV( a ) and 'Pa = kV( q ), where k is a constant, then we can show that 
K must be zero. (See Appendix 
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?V = 0, 
T aP c = 0, 

T ab c = 0, (4-10) 

and the equation of motion of supergravity, 10 ) 

T a p> = 0. (4-11) 

Using Eqs. (|4-1U|) and (|4-11|) . we have 

[V a , [V a , V 6 ]] = [V a , R ab cd O cd + T a ^V 7 ] 

= {V a R ab cd )O cd - R ab ac V c + (V a T a ^)V 7 + T ab -<R% cd O cd (4-l2) 

Furthermore, as shown in Appendix iBl from Eqs. (|4-1U|) . (|4-11|) and the Bianchi 
identities, we can show the relations 



Rab ac = 0, (4-13) 

V a T ab ~< = 0, (4-14) 

V a R ab cd = 0, (4-15) 

T ab -<R% cd = 0, (4-16) 



which indicate that the right-hand side of Eq. (|4-12|) is zero. Hence, the ansatz (|4-9|) 
with the constraints (|4-10j) and ()4-lljl satisfies Eqs. (|4-7j) and (|4-8j) . 

Some comments are in order here. First, although we showed that Eqs. (|4-10|) 
and (|4-11|) are sufficient conditions for Eqs. Q4-7|) and (|4-8|) under the assumption 
(14-91) . it is unclear whether or not they are necessary. Secondly, it is desirable to 
determine whether the torsion constraints in the ten-dimensional superspace 14 ) are 
compatible with the equation of motion of the supermatrix model. Thirdly, because 
we have set & a to zero as an ansatz, the action (j4-6|) allows the same classical 
solutions. 

4.3. Superdiffeomorphism and local supersymmetry 

The supermatrix model possesses the superunitary symmetry 

SAa = i[A,Aa], 5& a = i[A,V a ], (4-17) 

where A is an even Hermitian supermatrix. If we expand A a and !F Q as*-* 

A = iV (a) + a [a) (z,g) +ia ia) ( - B) {z,g)V {B) + ia {a) bc (z, g)O bc 

+i 2 a {a) {B){C) (z, g)V {B) V {C ) + i 2 a (a) ^ cd (z, g)V [B) O cd + ■ • • , 

& a = 4>( a )(z, 9) + i^(a) {B) (z, g)V {B) + iii (a) bc (z, g)O bc 

+i\a) mC \z,g)V {B) V (c) + i 2 ^ (a) ^ cd (z,g)V (B) O cd + • • • , 

V(A) = £(A) M Vm, (4-18) 

As in W2.'2\ because A a and Pa are Hermitian, we should introduce the commutator and the 
anticommutator. Here we omit them for simplicity. 
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and take A as 

A = \ {A< a >(*), V (a) } + \ [\ {a Hz), V (a) ] , (4-19) 
then Eq. (|4-17f) becomes a superdiffeomorphism on the fields that appear in Eq. 

gHEJ). 

The part of the superdiffeomorphism generated by 



A= % - 
2 L 



A (Q) (x),V (Q) J (4-20) 
gives the local supersymmetry. For example, in the Wess-Zumino gauge 8 )' 9 ) 



p A\ _ ( e m a (rr) 2 Xm a ( 3; ) i (a oi i 

^ |*=o " 1, J ' (4 21) 

where e m a is the vielbein and Xm a is the gravitino, Eq. ()4-20|) gives 

5e m a (x) = 2i\(x)-y a Xn(x), Sx m a (x) = 2V m X a (x) + (auxiliary fields)(4-22) 

§5. Conclusions and discussion 

In Ref. 1) we showed that the covariant derivative on a d-dimensional curved 
space can be described by a set of d matrices. Based on this, we introduced a 
new interpretation of the matrix model, in which we regard matrices as covariant 
derivatives on a curved space rather than coordinates. With this interpretation, 
the symmetries of general relativity are included in the unitary symmetry of the 
matrix model, and the path integral contains a summation over all curved spaces. 
In this paper, we have extended this formalism to include supergravity. If we pro- 
mote matrices to supermatrices, we can express supercovariant derivatives as special 
configurations of the supermatrices. Then, the local supersymmetry is included in 
the superunitary symmetry. Furthermore, we showed that if we consider the action 
S = —Str [AajA-b] 2 + • • • , the matrix equation of motion follows from the Einstein 
and the Rarita-Schwinger equations in the four-dimensional case. 

There remain several problems. First, it is desirable to clarify the relationship 
between our new interpretation and the original one, in which matrices represent 
coordinates. Although gravity is realized manifestly in the new interpretation, its 
relationship to string theory is rather obscure. Contrastingly, the original interpre- 
tation is directly connected to string theory, because the action can be regarded 
as the Green-Schwarz action of a type-IIB superstring. Furthermore, we have not 
understood the meaning of the global M = 2 supersymmetry of IIB matrix model 
in the new interpretation. Because we have realized supersymmetry as superunitary 
symmetry, it seems that global J\f = 2 supersymmetry need not exist. In the original 
IIB matrix model, however, it prevents the eigenvalues from collapsing to a point 
and allows a reasonable interpretation of spacetime. 3 ) We believe that in the new 
interpretation it plays a similar role and allows A a to fluctuate around a covariant 
derivative on some manifold. In order to connect the two different interpretations, 
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it may be helpful to consider noncommutative geometry. This is because there is 
no essential difference between coordinates and derivatives in noncommutative ge- 
ometry. 13 ) If we can understand the relationship more closely, the role of the global 
N = 2 supersymmetry may become clear. Furthermore, this would be useful for the 
construction of curved noncommutative manifolds. 

Secondly, it remains to investigate whether ten-dimensional supergravity can 
be derived from the equation of motion of our supermatrix model. In Refs. 14), 
the superspace formulation of ten-dimensional J\f = 1 supergravity coupled to super 
Yang-Mills was studied, and the torsion constraints were presented. The calculation 
in the ten-dimensional case is more complicated, due to the presence of various fields, 
including the dilaton, antisymmetric tensor and gauge field. At present, it is unclear 
whether the torsion constraints presented in Refs. 14) fit our supermatrix model. If 
they do not, we may need to modify the form of the action or to find other torsion 
constraints. 

Thirdly, the supermatrix model we have proposed might not be well-defined, 
because the action is not bounded from below, due to the minus sign appearing 
in the supertrace. The possibility that supermatrix models are formally equivalent 
to the corresponding ordinary matrix models is discussed in Refs. 15) and 16). *) 
If a similar mechanism holds in our supermatrix model, we may be able to embed 
supergravity in the original IIB matrix model. 
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Appendix A 

Supermatrices 

In this appendix we define a supermatrix and summarize its properties. 
Consider a supervector space V . We use the standard basis in which supervectors 
have the form 

v = ( Vl ) We , (A-l) 

V V 2 J }N ' y ' 

where v\ and V2 are bosonic and fermionic column vectors, respectively. We call v 
an even supervector. Multiplying v by a Grassmann-odd number, we obtain an odd 
supervector. 



*' As the action of the original IIB matrix model is formally the same as the effective action of 
D-instantons, our supermatrix model (14-211 is identical to the effective action of D-instantons and 
ghost D-instantons. 16 ' 
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In the standard basis defined above, an even supermatrix B and an odd super- 
matrix T can be written as 

N e No N e No 

Bi Fx \ }N e ( F[ B[ \ } N , 



s =u b.) n ) >»:■ (A - 2) 

Here I? and B' are bosonic matrices and F and -F' are fermionic matrices. An even 
supermatrix represents a linear transformation of the supervectorspace V, which 
maps an even (resp., odd) supervector to an even (resp., odd) supervector. An odd 
supermatrix also represents a linear transformation of V, but it maps an even (resp., 
odd) supervector to an odd (resp., even) supervector. 
A scalar supermatrix is defined by 

N e No N E No 

1 J }N ' V -1 J }N ' y J 

where c and e are Grassmann-even and odd numbers. Then, c commutes with both 
even and odd supermatrices, and e commutes (resp., anticommutes) with an even 
(resp., odd) supermatrix. 

The Hermitian conjugate is defined by 

--(sa)- *={%%)■ 

A supermatrix S is Hermitian if = S. Then, Hermitian matrices can be written 
as 

»=(3 bJ. ^=(3 I)' 

where Bi and F are bosonic and fermionic Hermitian matrices, and B and F are 
bosonic and fermionic complex matrices. 
The supertrace is defined by 



Str B = Str\^ ^ g ^ 1 = tr Si - tr £ 2) (A-6) 

5ir J" = Sir ^ ^} J = tr F{ + tr i^. (A-7) 

It satisfies the cyclicity relation 

Sir (5i5 2 ) = (-) SlS2 5tr («S 2 «Si) , (A-8) 

where s, is (resp., 1) if Sj is even (resp., odd). It also satisfies the following property 
for scalar supermatrices c and e: 

Str (pS) = c-Str S, Str (eS) = e ■ Str S. (A-9) 

Finally, the superdeterminant is defined for even supermatrices by 

sdetB = exp (Str log B) . (A- 10) 
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Appendix B 

Derivation of Eqs. \4-l°J[) , \4-14\) , \4'15\) and \4- lb\ l from Eqs. \4'1U[) and 

EH — 

In this appendix we show that if the torsion constraints and equations of motion 
of the four-dimensional supergravity in superspace, Eqs. (|4-1U|) and (|4-11|) . are sat- 
isfied, then Eqs. ()4-13|) . (|4-14() and ()4-15j) . and hence the equations of motion of the 
matrix model, are satisfied. The argument is almost parallel to that given in Ref. 
8). 

B.l. Bianchi identities 

The supertorsion and the super Riemann tensor satisfy the following Bianchi 
identities which arise from the super Jacobi identities. If we use Eqs. (|4-l(Jjl and 
()4-llj) to simplify the Bianchi identities, 
[V a , [V 6) V c ]] + (cyclic) = gives 

R[abc] = 0, (B-l) 
V [a R bc] ij +T [ai ?R c] v = , (B-2) 

[V a , [V 6 , V 7 ]] + (cyclic) = gives 

-R^J + Ra^ - 2iT ab s ( 7 l ) 75 = 0, (B-3) 

(7ii)/ - V 7 T a6 5 = 0, (B-4) 
and [V a , {V^, V 7 }] + {V^, [V 7 , V a ]} - {V 7 , [V a , V^]} = gives 

Rfr ij = 0, (B-5) 

2i (j k ) R ak ij - VpR aj ij - V 7 iV j = 0, (B-6) 

2i (j k )^T ak s - Y -R a ji ( 7 ii)/ - \R a p lj (7«)/ = 0, (B-7) 

where we have used Eq. ()B-5j) to simplify Eq. ()B-6(I . Here we have presented only 
the identities that we need in the following. 

B.2. Derivation of Eqs. WW , WW and 

B.2.1. Proof of R ab cd = R cd ab 

From the fact that R a b cd is antisymmetric under the exchanges a <-> b and c <-> d 
and Eq. HB-lj) . it follows that 

Rabcd — Rbcad Rcabd 
= —Rcbda + Rcadb 
— Rdcba ~i~ Rbdca Radcb Rdcab 
= ZRcdab + Rdbac + Radbc 
= 2R C( iab Rbadc 

= ZRcdab - Rabcd- (B-8) 
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Hence, we have 



ttab — -K ab- 



B.2.2. Proof of the Rarita-Schwinger equation 
Multiplying Eq. (|B-7|) by j(7fc) 7/3 , we have 

1 



2iT ab ° = -Ra^frbJij) 



l3S 



= \ R aP 3 (Vbijj ~ Vbjli + Ibijf 6 
= \Ra^{liY & + \Ra^(lbi J Y 5 

= \ V(7,) 5/3 " | (B-10) 
Then, multiplying this by (7 fe ) 7 5, we obtain 

2z^ b rsT ab S = \Rap bj (l bj yP - \Ra^{l b lU^ 

= \^ r'i'ur ! - ^ V'(7^) 7/3 

= 0. (B-ll) 

On the other hand, by contracting 7 and 8 in (|B-7|) . we have 



2i^ b r S T ab s = ^Ra^hv)"" 3 - (B-12) 

Therefore, we find 

(l a ) aP Tj = 0, R a p ij {lijV P = 0. (B-13) 
The former is the Rarita-Schwinger equation written in terms of superspace. 8 ) 
B.2.3. Proof of the Einstein equation (|4-13|) 

By contracting a and i in Eq. (|B-3|) and using the second equation in Eq. ()B-13|) , 
we find that 

R ai ac = 0. (B-14) 



Next, contracting a and i in Eq. ()B-6|) and substituting the above equation, we 
obtain 

(7 k )(3 7 Rak aJ = 0, (B-15) 

which is equivalent to Eq. (|4-13|) : 

Rab ac = 0. (B-16) 

This is the Einstein equation without a cosmological constant. Note that the su- 
per Ricci tensor R a b ac contains the contribution from the energy-momentum of the 
gravitino. 
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B.2.4. Proof of Eq. ifTOl) 

By Eqs. (fTTU)) and (|B-5I) . we have 

{V Q ,V^} = 2i( 7 a C- 1 )^V a , (B-17) 
V a = -^(C 7 ar /3 {V a ,V /3 } 

o 

= ~{C-y a )^V a Vp. (B-18) 
Then, multiplying Eq. (|B-4|) by ( 7 a )p 7 , we have 

1 



= 0. (B-19) 

Here we have used the relation R a \ ba q = 0, which follows from Eqs. (|B-9|) and (|B-1|) . 
and Eq. I|B-16|1 . Next, multiplying Eq. (|B-19|) by V p and using Eq. (|B-18|) . we 
obtain Eq. (|4-14|) : 

V a T afe <5 = 0. (B-20) 

B.2.5. Proof of Eqs. (fIT5j) and (I^Tfil) 

From Eqs. (lB~T3|) . (IB~3l) and (lB~7l) . it follows that 8 )' 9 ) 

iU cd = 2i( lb C- 1 UT cd(3 . (B-21) 
Then, because C 7 a is symmetric, we have 

T bc ^R a ^ = T l ^R a lbc . (B-22) 
Combining this with Eq. ()B-9|) . we can rewrite Eq. ()B-2|) as 

V [a Rij bc] + TifR^J^ = 0. (B-23) 

The second term on the right-hand side is zero, because substituting Eq. (|B-21|) into 
Eq. (|B-3|) we have R^^ = 0. Therefore, we obtain 

V [a Rij bc] = 0, (B-24) 

and contracting a and i, we have 

V a R ab cd + V c R ab da + V d R ab ac = 0. (B-25) 
Because the second and third terms are zero, by Eq. (|B- 16|) . we obtain Eq. (|4-15|) : 

V a R ab cd = . (B-26) 

From Eq. (lB~2ll) we have 

= T ab ~< ■ (-2i)( 7 a C- 1 ) 7(5 r i ^ 

= 2*2** ■ (^C- 1 )^^, (B-27) 
but this is zero, by Eq. (|B-13|) . Therefore, we have shown Eq. (|4-16|1 . 
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Appendix C 

Proof of k = 

In the footnote of ^4,21 we claimed that if we stipulate A a = ^ (a) > = «V( Q ) 
and the standard off-shell torsion constraints ()4-lU(l as an ansatz, then the equations 
of motion (|4-7|) and 1|4-8|) of the matrix model force k to be zero. In this appendix, 
we give an outline of its proof. 

By combining Eq. (|4T0j) with the Bianchi identities, the remaining degrees of 
freedoms can be expressed in terms of three superfields: the complex scalar R, the 
real vector G a , and the spin | field W( a p y y 9 ^' 12 ^ 

If K 7^ 0, Eq. (H33> reduces to 

(7 a )/ (R a p cd O cd + 7VV 7 ) = 0, (C-l) 

which is equivalent to 

{l a )a P Ra P Cd = 0, (C-2) 

(-/V T af ? = 0. (C-3) 

In terms of G, R and W, the latter becomes *) 

G a = R = 0. (C-4) 

This is the equation of motion of supergravity, 10 ^ and it is equivalent to Eq. (|4-llj) . 
As explained in the previous subsection, this implies Eq. (|B-16j) : that is, the cosmo- 
logical constant is zero. (Furthermore, Eq. (|(J-2|) gives VF( Q/ g 7 ) = 0, which with Eq. 
(|C-3|) allows only a flat spacetime.) On the other hand, substituting Eqs. (|4-9|) and 
(|4-1U(I into Eq. (|4-7|) . the coefficient of V c becomes 

i {R ab ac + 8kV) = 0. (C-5) 

Hence we find that the cosmological constant is given by A = -32k 2 . This is a 
contradiction, and thus k must be zero. 
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